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Cosmological models involving a bounce from a contracting to an expanding uni- 
verse can address the standard cosmological puzzles and generate "primordial" den- 
sity perturbations without the need for inflation. Some such models, in particular 
■ the ekpyrotic and cyclic models that we focus on, fit rather naturally into string the- 

ory. We discuss a number of topics related to these models: the reasoning that leads 

& ■ 

to the ekpyrotic phase, the predictions for upcoming observations, the differences 
between singular and non-singular models of the bounce as well as the predictive 
and explanatory power offered by these models. 
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I. INTRODUCTION 

Our universe is special in that it evolved from an a priori highly unlikely state, in the sense 
of Boltzmann: the very early universe was exceptionally flat, homogeneous and isotropic 
(hence it had very low entropy), with approximately scale- invariant density perturbations of 
magnitude Q ~ 10~ 5 added on top. Now it may be that string theory, once better under- 
stood, leads to a unique initial state with these characteristics, but there is currently little 

O ■ 

evidence pointing in that direction. Therefore we would like to find a dynamical explanation 
for these "initial" conditions. 
. The most popular attempt at formulating such a model is inflationary cosmology [1- 

3]. In recent years, many attempts have been made at obtaining inflationary dynamics in 



string theory (see for example the contributions by McAllister and Burgess, Quevedo and 
Mazumdar to this focus issue). All models to date require rather elaborate constructions 
(with multiple branes of diverse dimensionalities, orientifold planes and warped geometries 
induced by p-form fluxes all carefully balanced) together with their own set of initial con- 
ditions. Thus, currently, these models simply report the question of initial conditions to an 
earlier epoch, perhaps alleviating it, but not solving it. What these models do achieve is that 
they provide a mechanism with a microphysical basis for generating the primordial density 
perturbations, and this is their main merit. All models to date still have to be fine-tuned 
in order for the characteristics of the generated perturbations, such as their amplitude and 
spectrum, to agree with observations (see also the contribution by Mulryne and Ward [4]). 
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This unfortunately reduces the explanatory power of inflationary models. A recourse that is 
sometimes advocated is to invoke eternal inflation [5] (which creates an infinite number of 
universes of all possible configurations - hence it will also create universes with the elaborate 
characteristics necessary for ordinary inflation to occur) together with anthropic arguments 
[6] (which "justify" disregarding all universes with characteristics substantially different from 
ours). Such an approach may however turn out to be fallacious: eternal inflation leads to 
infinities, which must be regulated. However, all results turn out to depend on the regulator, 
i.e. eternal inflation does not make regulator-independent predictions (this is the so-called 
measure problem of inflationary cosmology, see the contributions by Freivogel [7] and Kleban 
to this focus issue). Hence one may suspect that eternal inflation is not a consistent theory 
1 , and we will adopt this point of view in this article. In this respect, it certainly seems 
worthwhile to examine alternative dynamical models for the early universe. In this article, 
we will discuss models that are based on a cosmic bounce, i.e. models in which the current 
expanding phase of the universe was preceded by a contracting phase. We will see that such 
models, in particular the ekpyrotic [10] and cyclic [11] models of the universe, may provide 
a framework that can explain the "initial" conditions of the universe dynamically, without 
recourse to anthropic arguments. Moreover, these models make rather distinctive predictions 
for upcoming cosmological observations, which we will discuss. 

The original ekpyrotic model was inspired by string theory [10] , in particular by Hofava- 
Witten theory [12, 13]. We start by reviewing a simplified version of the string theoretic 
setup in question (section II), and build up the logic that leads to a cyclic universe from 
there. In brief, the logic goes as follows: the big bang may be seen as a collision of branes in 
higher dimensions, but where the scale factor on the branes does not shrink to zero size, only 
the distance between the branes goes to zero. Hence the big bang may have been a much 
milder, and potentially tractable, event. If so, the universe would have been in a contracting 
phase previously. However, a contracting universe is unstable to the buildup of anisotropies 
and the onset of chaotic mixmaster behavior, unless there is a matter component in the 
universe with ultra-stiff equation of state. This matter component is "ekpyrotic" matter 
(section III), and it can be modeled as a scalar field with a steep and negative potential (this 
scalar may be the scalar parameterizing the distance between the branes, its potential would 
then describe an attractive inter- brane force). Surprisingly, this ekpyrotic phase not only 
manages to flatten and isotropize the universe, but it also amplifies quantum fluctuations 
which may obtain a scale-invariant spectrum via the entropic mechanism reviewed below 



Albrecht [8] and Banks [9] have suggested that holography may prevent eternal inflation from occurring. 
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(section IV). As the potential is steep, the scalar is strongly self-coupled, and the resulting 
perturbations are distinguished by large non-gaussianities of the local form. However, as 
for inflation, for the ekpyrotic phase to be successful, the universe must again have had the 
right initial conditions, so it looks like we might not have made any progress at all. This is 
where the idea of a cyclic universe comes in (section VI): it turns out that by joining the 
ekpyrotic phase to the future of the current dark energy phase, it becomes manifest that the 
current dark energy phase sets up precisely the right conditions for a subsequent ekpyrotic 
phase to be successful. This then leads to the idea of a cyclic universe, in which the most 
recent ekpyrotic phase (during which our perturbations were generated) was also preceded 
by a dark energy phase etc. Towards the end of this article, we will present different models 
for the all-important bounce phase (section V) and discuss the reasons why such a cyclic 
framework may be powerfully predictive (section VI). 



Consider pure gravity in 5 dimensions. Motivated by Hofava-Witten theory [12, 13], we 
want to look for a solution where one spatial dimension consists of a line segment. A line 
segment can also be described as the quotient, or orbifold, of a circle by a Z 2 reflection 
symmetry across a diagonal of the circle. Under the action of the Z 2 reflection symmetry, 
the two endpoints of the line segment are fixed points - from the 5-dimensional point of view 
these define the location of two orbifold planes, which are (3 + l)-dimensional boundaries of 
spacetime. We can write the metric as 



where — < y < i/q denotes the orbifold coordinate and is the radion field parameterizing 
the size of the orbifold. The ^-dependent prefactor in front of the 4-dimensional metric 
&s\ ensures that after dimensional reduction we are left with the canonically normalized 



II. COLLIDING BRANES AND THE BIG BANG 




(1) 



Lagrangian 2 




(2) 



Thus, the equations of motion become 



3H 2 = -ft 
2 



-H 



(j) + 3H(f) = 0. 



(3) 



2 We use reduced Planck units h = c = 1 and 8nG = M p ? = 1 . 
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They imply that a oc e' 



and are solved by 



a = a (-t) 1 / 3 




ln(-t) + 0o, 



(4) 



for some integration constants a , 0o- This means that we can now plug the solution (4) into 
the original 5-dimensional metric, to find (up to obvious re-scalings) 



where we have changed the time coordinate to T oc (— t) 2 ^ . This model spacetime, sometimes 
called the compactified Milne spacetime, describes two (empty) orbifold planes approaching 
each other, colliding, and receding away from each other again. Under the further change 
of coordinates u = Tcoshy, v = Tsinhy, the metric becomes simply that of 5-dimensional 
Minkowski space 



which shows that this spacetime is flat, except at the moment of collision (T = 0), see 
also figure 1. From the 4-dimensional point of view, the moment of collision corresponds to 
a{t) = 0, i.e. it corresponds to the big bang. However, as seen from the full 5-dimensional 
theory, the singularity corresponds to only the orbifold dimension shrinking to zero size, and 
not the other spatial dimensions. Indeed, as can be seen from (1), the brane scale factors 
are really given by e~^^a, and they are constant and non-zero at the collision. The density 
of matter (which is stuck on the branes and subdominant regarding the dynamics of the 
collision) as well as the temperature remain finite at the collision. 

This simple example shows that, in a higher- dimensional context, the big bang singularity 
may be much milder than in 4-dimensional general relativity, and may be tractable. We will 
discuss models of cosmic bounces (both singular, like this one, and non-singular) in more 
detail in section V. For now, let us just go with the idea that the big bang may not have 
been the beginning of the universe, and that there may have been a previous phase during 
which the universe contracted. This immediately leads to a potential problem, namely the 
wild growth of anisotropies, combined with a chaotic mixmaster big crunch, that is expected 
to occur in a contracting universe. In such a case, we evidently wouldn't expect a flat and 
isotropic universe to re-emerge from the crunch. One way to avoid this would be for the 
universe to be exceptionally flat at the start of the contracting phase, but this would amount 
once more to imposing special initial conditions - the very thing we want to explain. But 
there is another possibility, described in [14], which involves the addition of a new matter 
component. We will discuss it next. 



dsl = (-t)- 2 / 3 [-dt 2 + (-t) 2 / 3 dxj} + (-tf 3 dy 2 
= - d T 2 + T 2 dy 2 + dx 2 3 , 



(5) 
(6) 



dsj = ~du 2 + dv 2 + dx 2 3 , 



(7) 
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FIG. 1: The compactified Milne universe (the shaded region) is a useful model for a big bang caused 
by the collision of two boundary branes. This figure illustrates the embedding in 5-dimensional 
Minkowski space, the 3 ordinary spatial dimensions being suppressed and with the dashed lines 
delineating the light cone at the moment of the collision. 

III. THE EKPYROTIC PHASE 

For concreteness, consider a Friedmann-Robertson- Walker (FRW) metric 

\ 1 — nr z J 

where a(t) denotes the scale factor of the universe and k = —1, 0, 1 for an open, flat or closed 
universe respectively. If the universe is filled with a number of fluids interacting only via 
gravity and with energy densities pi and constant equations of state tOj, then the equations 
of continuity 

Pi + 3-( Pi +Pi) = (9) 
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(where dots denote derivatives with respect to time t) imply that they will evolve according 
to 

Pi cx a~ 3 ( 1+Wi \ (10) 

The Einstein equations lead to the Friedmann equation, which involves the Hubble parameter 
H = a/a: 

H 2 = U^ + ^f + ^f + ^f + ... + ^^). (11) 

3 V a 2 a 3 a 4 a 6 a 3 ( 1+ ^) / v ; 

The p^s are constants giving the energy densities at scale factor a = 1 of the various 
constituents of the universe: non-relativistic matter (subscript to), radiation (r) and the 
energy density associated with anisotropies in the curvature of the universe (a). There is 
also a term due to the mean curvature of space, and we assume the presence of a scalar field 
0. 

In a contracting universe, as a — > 0, the term scaling with the largest negative power of 
a will eventually come to dominate. And this will be the anisotropy term oc a -6 , unless 
the scalar field has an equation of state w$ > 1. In that case, the scalar field will come to 
dominate and suppress the anisotropies, as well as the mean curvature, matter and radiation. 
Such a phase is called an ekpyrotic phase. Since the equation of state of a scalar field is given 
by 

P y> 2 - v{4>) 

w d> — ~ — "T^ > (12) 

* P |0 2 + \/(0) 1 ^ 

we can see that w ( p> 1 can be achieved by having a negative potential V((f>). A simple, but 
relevant example is the negative exponential 

V(4>) = -Voe-* (13) 

where V and c are constants. As mentioned previously, in string theory such a scalar could 
be the radion, i. e. the scalar field describing the distance between the two boundary branes, 
and in that case a negative exponential potential can be generated by instanton effects, see 
e.g. [15, 16]. Once the scalar comes to dominate the dynamics, the equations of motion 
reduce to 

3H 2 = \$ - Voe- * <j> + - cV^T^ = 0, (14) 
and they are solved by 

a = (-tf/\ 0=y|ln(-v^t), € =j. (15) 



The equation of state is 



2e c 2 

W * = J- 1 = J- 1 > ( 16 ) 
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so we need c > y/6 for ekpyrosis. In realistic models (realistic in the sense of agreeing 
with observations), e typically turns out to be of 0(100). Note that the same parameter, 
in inflationary models, is the "slow-roll" parameter and is typically of 0(1/100). Here e is 
rather a "fast-roll" parameter. Note that when e is large, the potential is steep and the 
universe is contracting very slowly. 

IV. PREDICTIONS FOR OBSERVATIONS 

As described above, the ekpyrotic phase suppresses the mean curvature and anisotropies, 
and so, if this phase lasts long enough (typically a fraction of a second is already enough, see 
e.g. [17]), the flatness problem is solved. Incidentally, the horizon problem is automatically 
solved too, as, during the contracting phase, there is ample time for the different regions we 
observe in the cosmic microwave background (CMB) to have been in causal contact. However, 
the ekpyrotic phase not only resolves the standard cosmological puzzles, it also generates 
density perturbations by amplifying quantum fluctuations [18], much like inflation, but with 
some crucial differences in the details. Here we simply summarize the results - the detailed 
derivations were recently reviewed in [19]. 

During the ekpyrotic phase, the universe is contracting very slowly, but the horizon 
1/ H cat — 7-0 is shrinking fast. Hence, similarly to what happens during inflation, quantum 
fluctuations with progressively smaller wavelengths exit the horizon and turn into classical 
density perturbations. Now it turns out that if there is a single scalar field, with a steep 
and negative potential, driving the ekpyrotic phase, then the resulting curvature fluctuations 
have a blue spectrum (with spectral index n s « 3) [18, 20-22] 3 . On large scales, these per- 
turbations are completely subdominant (although they could make a re-appearance during 
the bounce phase, as discussed further on) and hence these cannot give rise to the temper- 
ature fluctuations we observe in the CMB. But if there is more than one ekpyrotic scalar 
present, scale-invariant curvature perturbations can be generated via the entropic mecha- 
nism [25-27]. In fact, from the point of view of string theory, it is very natural to expect 
more than one scalar. For example, in the colliding branes setup described above, there was 
one scalar field describing the distance between the two orbifold planes along a fifth dimen- 
sion. However, in the full string/M-theory setting, we would expect there to be 6 additional 
internal dimensions, for example in the shape of a Calabi-Yau manifold [28]. The volume of 

3 Right at the onset of the ekpyrotic phase, a small range of scale- invariant curvature perturbations are 
produced by the adiabatic mechanism [23, 24]. However, it has so far not been possible to incorporate this 
mechanism into a full ekpyrotic or cyclic model of the universe. 
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the Calabi-Yau manifold is then described by an additional scalar, and many more scalars 
are required to describe all the shape deformations of the internal space. Since the inclusion 
of one additional scalar already entails the needed qualitative differences to the single-field 
case, we will focus on having two scalars for simplicity, with action 

S = J ^[R- l -{d4>,f - l -{d<t> 2 f -V^^)]. (17) 

Let us assume also that both fields get an ekpyrotic-type potential, 

Vfa, <t> 2 ) = -V ie - C ^ - V 2 e~ c ^\ (18) 

Then it is much more natural to discuss the dynamics in terms of the new variables a and 
s pointing transverse and perpendicular to the field velocity respectively [29, 30]; they are 
defined, up to unimportant additive constants which we will fix below, via 

0101 + 0202 0102 - 0201 n n , 

a = : , S = : , (19) 

a a 

with cr = (0f + 4>\) l l 2 . For later use, we also define the angle 9 of the trajectory in field space, 
via [31] 

cos^ = — , sin# = — . (20) 
a cr 

In terms of these new variables, the potential can be re-expressed as 

Vek = -V e^[l + es> + |e 3 /V + ^eV + •••], (21) 

where for exact exponentials of the form (18), one has k 3 = 2\^2(cf — c|)/|cic 2 | and 
k 4 = 4(cf + cf) I (cfcKcl + c|)). However, in the absence of a microphysical derivation of 
the potential, we will simply take k 3 ,K4 ~ 0{\) and express all results in terms of «3,«4. 
The ekpyrotic scaling solution becomes 

a(t) = (-t) 1/t a = In (-y/^tj s = 0, (22) 

with the angle 9 being constant. Thus we can see that, up to a sign determined by convention, 
the field a now plays the role of the single scalar we discussed earlier, and hence this is the 
field that picks up a blue spectrum of curvature perturbations. The transverse field s, on 
the other hand, feels a different potential. Because of the overall minus sign in (21), we can 
see that the potential for s is actually unstable [27, 32], and the motion during the ekpyrotic 
phase is along a ridge in the potential. This has profound consequences, which we will discuss 
in section VI. In the present context, the instability implies that the s field fluctuations Ss, 
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which are entropy (or isocurvature) fluctuations, get amplified. More detailed considerations 
show that they obtain an amplitude Q s and a spectral index n s given by [27] 

Q s ~ {eV mm ) 1 ' 2 (23) 
n,-l = H_f£, (24) 

where V m i n denotes the value of the ekpyrotic potential at the end of the ekpyrotic phase. 
The fast-roll parameter e is now allowed to vary slowly and diV = din a. Hence, for e of 
(9(10 2 ) (or more), the spectrum is close to scale- invariant, with the first term on the right- 
hand side tending to make the spectrum blue and the second term tending to make it red 4 . 
A simple estimate, based on [33], of the natural range of n s gives the range 0.97 <n s < 1.02, 
in good agreement with current data. 

In this way the ekpyrotic phase generates nearly scale-invariant entropy perturbations. 
But what we are really interested in are the curvature perturbations (, which are the per- 
turbations that give rise to the temperature fluctuations in the CMB. In that respect it is 
useful to examine the evolution equation for the (linear) curvature perturbations on large 
scales, given by [31] 

C = ~^y^ s - (25) 

As soon as the background trajectory bends (6 ^ 0), the entropy perturbations source 
the curvature perturbations, implying that on large scales the curvature perturbations thus 
created obtain the same spectrum as the entropy fluctuations (there is no wavelength de- 
pendence in (25)), with an amplitude of [19] 

Ql - (26) 

We will discuss the amplitude in much more detail in section VI. In the most concrete models 
of the entropic mechanism studied in the literature (reviewed in section V), the indispensable 
bending of the trajectory occurs automatically. These models show that the bending occurs 
after the ekpyrotic phase has come to an end, in the phase leading up to the bounce, when the 
ekpyrotic potential has become unimportant and the dynamics is dominated by the kinetic 
energy of the a field 5 . 

4 Note that during the ekpyrotic phase, ./V is decreasing, and so must e in order for the ekpyrotic phase to 
come to an end eventually. 

5 It was also previously envisaged that the trajectory could bend due to the instability of the potential 
[29, 34]. This case is now disfavored, as it leads to predictions for non-gaussianities in disagreement with 
observations at the 4ct level, see [19, 35]. 
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The curvature perturbations created in this way are indistinguishable from those created 
by a phase of inflation, for example. However, at higher orders, the ekpyrotic perturbations 
are characterized by sizeable and specific non-gaussian signatures [19, 35-37]. There is 
a simple reason for this: the potential is necessarily steep during ekpyrosis, so that the 
scalar fields are necessarily significantly self-coupled, leading to non-gaussian statistics for 
the perturbations. It turns out that these non-gaussian corrections are of the so-called "local" 
form. If the full, non-linear curvature perturbations ( are re-written as a linear, gaussian 
piece Q, plus correction terms, 

C = C 9 + \InlQ + ^9nlC 3 9 , (27) 

then the predictions for ekpyrotic curvature perturbations generated via the entropic mech- 
anism are [38-41] 

Inl = ^3V^+5 (28) 

9nl = (|« 4 + ^-40)e. (29) 

With k 3 ,k 4 expected to be of 0(1), and e typically of C(10 2 ) (e needs to be above about 
50 in order for the power spectrum to be in agreement with observations), we would thus 
expect the parameter f NL to be of C(±10), with the sign being typically determined by the 
sign of k 3 , and g^L of 0(— 1000) and negative in sign. The current observational bounds 
are that f NL = 38 ± 21, where the errors are quoted at one standard deviation [42], while 
currently no strong constraints exist on g^L. It will therefore be extremely interesting to 
see the improved data likely to be gathered by the PLANCK satellite and ongoing and 
future large scale surveys, especially in light of the fact that the ekpyrotic predictions are 
significantly different than those of simple inflationary models. 

Before proceeding to discuss the bounce in more detail, we should also discuss the predic- 
tions for primordial gravitational waves [10, 43, 44]. In fact, a quick heuristic argument is 
sufficient to get at the main result: at the linear level, the generation of gravitational waves 
depends entirely on the behavior of the scale factor. And during the ekpyrotic phase, the scale 
factor shrinks very slowly, so that, to a first approximation, the background is Minkowski 
space. This immediately implies that no significant large-scale gravitational waves are gen- 
erated during ekpyrosis (otherwise gravitational waves would also be produced all around us 
all the time), in stark contrast to large-field inflationary models. It turns out that, at second 
order in perturbation theory, small-amplitude gravitational waves are sourced by the scalar 
density perturbations, and these are the dominant gravitational wave signal expected from 
an ekpyrotic phase [45], though they are unlikely to be observable in the near future. Hence 
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a detection of large-amplitude primordial gravitational waves would rather conclusively rule 
out the models described here. 

These predictions are interesting, because they are quite specific: reasonably large /atl, 
reasonably large and negative g^L plus an absence of primordial gravitational waves. It is 
clear that inflationary models can also be constructed with exactly these predictions, hence, 
even if future observations remain in agreement with the predictions, no clear-cut conclusions 
can be drawn. However, inflationary models with exactly these predictions would have 
to be rather contrived and fine-tuned, and would have to be adjusted a posteriori, which 
would render them disfavored in the author's opinion (and according to standard results in 
probability theory). 

V. COSMIC BOUNCES: SINGULAR OR NON-SINGULAR? 

In cosmological models involving a reversal from contraction to expansion, the bounce 
phase is the most crucial ingredient. This is the moment that gets identified with the "big 
bang", and it is the time when quantum gravity effects are likely to play a dominant role. 
Since we currently do not possess an experimentally supported theory of quantum gravity, 
all current models are necessarily rather speculative. Nevertheless, using current candidate 
theories, in particular string theory, we can construct and analyze various models for a 
bounce, and explore their consequences. 

During the contracting phase, the Hubble rate is negative, while it is positive during the 
subsequent expanding phase of the universe. Hence, the bounce must allow the Hubble rate 
H to increase. If we describe the dominant energy component in the universe as a perfect 
fluid with energy density p and pressure p, then the Einstein equations imply that 

H = - 1 -(P + P). (30) 

Obtaining the required increase in the Hubble rate can then happen in only two ways: either 
p+p < 0, i.e. the null energy condition is violated, or the bounce must be classically singular 
with the Hubble rate instantaneously jumping from negative to positive values. Of course, 
especially in the second case, it becomes crucially important that the quantum corrections to 
the classical description be understood. We will discuss models of both types below. Before 
doing so, it may be useful to point out that the case of a singular bounce is not a limit of 
a classically non-singular bounce. In the non-singular case, the scale factor of the universe 
must slow down its contraction, come to a halt, then increase again. This is typically a very 
fine-tuned process. By contrast, in the singular case, there is typically no slowing down, as 
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will be evident in the example described below. Hence we should not expect these two types 
of models to lead to equivalent physical results. 



A. Classically Singular Bounces 

The most interesting example of a classically singular bounce is the collision of the orbifold 
branes in Horava-Witten theory/heterotic M-theory, a simplified version of which was already 
described in section II. Here we will elaborate on this process, and list the reasons for why 
this example is particularly encouraging. 

Horava-Witten theory describes the strong coupling limit of the 10-dimensional E s x E$ 
heterotic string theory [12, 13]. The value of the string coupling is determined by a field 
(the 4> field of section II) which gets identified with the size of a new spatial dimension. As 
the coupling constant grows, this eleventh dimension grows too, and in this newly opened 
space the theory becomes that of 11-dimensional supergravity. The endpoints of the new 
11th dimension consist of two orbifold fixed planes, which are endowed with an E 8 gauge 
theory each. If in addition six of the spatial dimensions are compactified on a Calabi-Yau 
manifold, then the theory, now usually called heterotic M-theory, contains (3 + l)-dimensional 
boundary branes with realistic particle physics on them [28, 46, 47]. Restricting our attention 
to gravity and the Calabi-Yau volume modulus Vcy only (these are the only fields that are 
relevant close to the bounce), the action is 



J5d 



R ~ \v^d m V CY d m V C Y - 6a 2 V£ 



+-L(-12a / V=9~Vc Y + 12a [ V=^V^\ , (31) 

ZK 5 I J4d,y=+1 J4d,y=-1 ) 

where a denotes the boundary brane tension. The static vacuum solution to this theory is 

ds 2 = h 2 /\y) [A 2 (-dr 2 + df 2 ) + B 2 dy 2 ] , 
Vcy = Bh % /\y), 

h(y) = 5ay + C, -1 < y < +1 (32) 

where A, B and C are integration constants. The y coordinate spans the orbifold with fixed 
points at y — ±1 (we have re-scaled y compared to section II), a negative-tension brane 
being located at y — — 1 and a positive-tension brane at y — +1. As we will justify below, 
we are interested in the situation in which the branes are moving slowly. In order to describe 
that case, we can simply promote the moduli A, B, C to functions of (conformal) time r. 
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Then, after the field redefinitions 

a 2 = A 2 Bh, (33) 

5 

D ^/V2 = B{Izf/\ (34) 

5 



02 



20 



j dC{Iz)- l ^{I_^ + lQI_ l I^ , \ (35) 



dy h» = - . — [(C + 5«)(" +1 ) - (C - 5«)(" +1 )], (36) 
5a:(n + 1) 

the action simplifies considerably, to [48] 

S mod = / [-3a' 2 + ia 2 (0; 2 + 0; 2 )], (37) 

where ' = This is simply gravity with scale factor a minimally coupled to two scalars, 
just like the theory we used to describe the entropic mechanism in section IV, except that 
we have written it our here on a FRW background and in terms of conformal time. Higher- 
dimensional quantities are then determined in terms of the fields a, 0i,02, for example the 
distance d between the branes and the values of the brane scale factors a± at the locations 
of the boundaries at y — ±1: 

d = L_ e <h./V2-y/3/24>2 | c 2 A /273<fe)3/2 ^ ^/iT^fe 1 3/2j ^ ^gg^ 

1/4 




cosh a/2/302 

a ± = (2a)^ae-^ 2 <J ) 7 , i/4 (39) 

^- sinh a/2730 2 J . 

The axis 02 = is special: positive values of 2 would lead to unphysical imaginary values 
for the brane scale factor a_, and are thus forbidden. One can show that in the presence 
of matter on the boundary branes, there arises an effective repulsive potential along the 02 
axis, which causes the field space trajectory to bend in its vicinity [48, 49]. This is of course 
just what is required for the conversion of entropy into curvature perturbations described in 
section IV, and the example here provides a microphysical basis for the entropic mechanism. 

The effective theory (37) contains two scalar fields, unlike the toy model of section II. 
However, by requiring that at the brane collision the brane scale factors do not shrink to 
zero or blow up to infinity, we are forced to look at solutions where close to the collision 
0i — v / 302 approaches a constant. This effectively reduces the number of scalars to one, when 
looking at the background evolution (for the perturbations, it remains of course important 
that the theory really contains two scalars). The theory (37) is then solved by 

a = \2y r\^, = \2y Q rf^\ = \4ay r\^\ (40) 
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where 2y has the interpretation of being the velocity of the branes at the collision, as (38) 
implies that d ~ 2y r for small r. Note that yo is the only parameter of the solution. 

Close to the collision, the brane tension a becomes irrelevant and these solutions approach 
the compactified Milne solution described in section II, with the speed at collision being given 
by 2y . Their most important feature is that in the approach to the collision the spacetime 
that they descibe becomes increasingly flat (of course, at the moment of collision there is still 
a 5-function curvature singularity). Hence, higher-derivative corrections induced by quantum 
loops become more and more negligible too [50]! This is similar in spirit to a conjecture due 
to Penrose, which states that quantum corrections are small when the Weyl curvature is 
small (as is the case here too) [51]. So, even though this solution breaks supersymmetry 6 , 
the quantum corrections can be expected to be mild. This hope is given further support by 
the fact that the distance between the orbifold branes also determines the string coupling 
constant, as seen from 10 dimensions. As the branes approach each other, this coupling 
decreases, and it is zero at the moment of collision. Thus, during the bounce, we expect 
interactions to be significantly suppressed. This can be made more quantitative by studying 
certain aspects of the brane collision semi-classically. In particular, membranes that wind 
around the compact dimension are expected to become light as the size of the eleventh 
dimension shrinks. Thus such states must be taken into account close to the collision. Turok 
et al. have shown that the equations of motion for winding M2-branes remain regular at 
all times, and that they propagate smoothly through the collision [53] - see also [54-56]. 
Instanton techniques then enable one to estimate the production of such light membranes 
at the collision. The result is that the density of produced membranes is finite and their 
gravitational backreaction is small, provided the collision speed yo remains non-relativistic. 
This last finding has interesting consequences regarding the amplitude of the cosmological 
perturbations that are produced in these models, as will be discussed in the final section 
below. 

Eventually, one would like to be able to treat the collision of the orbifold planes of heterotic 
M-theory fully quantum mechanically. Unfortunately, this seems to require that a full non- 
perturbative definition of M-theory be found first. This drawback notwithstanding, it is 
interesting that all indications found so far imply that this particular example of a singular 
bounce may well be viable. 

6 Even the exactly flat compactified Milne solution breaks supersymmetry, because the fermions are not 
invariant under the boost identifications that need to be performed in order for one dimension to be a line 
segment [52]. 
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B. Classically Non-Singular Bounces 



As we saw at the beginning of this section, in order for the universe to bounce in a 
non-singular way, the dominant energy component in the universe must violate the null 
energy condition (NEC). From a model-building point of view, this is difficult to achieve, as 
violations of the NEC are commonly associated with various instabilities in the theory (most 
notably the appearance of ghosts, i.e. fields with negative kinetic energy). Moreover, since 
the energy density scales as a~ 3( - 1+ "'* ) for matter with equation of state w, and since violations 
of the NEC correspond to w < — 1, we can see immediately that in a contracting universe 
such matter becomes quickly less and less relevant. In fact, as discussed in section III, at 
the end of the ekpyrotic phase all other types of matter have become vastly subdominant 
to the ekpyrotic scalar field. Hence, the only option seems to be that the ekpyrotic scalar 
itself must turn into the required NEC- violating matter [34]. At that point, our options are 
currently rather limited: there exist only two types of scalar field theories in which the NEC 
can be violated in a stable and controlled way: ghost condensates [57] and galileons [58-60]. 7 
We will mainly focus on the most studied case, that of a ghost condensate bounce, which 
has been used in "new ekpyrotic" models [34], see also [37, 62]. It works as follows: suppose 
that, as the kinetic energy for the scalar field becomes very small, its kinetic term becomes 
non-canonical. In fact, if we denote -^r^) 2 = AT, so that the ordinary kinetic term in 
the Lagrangian is simply A, then we imagine that, for small A, the Lagrangian (which we 
denote by £ = -P(A)) is of the form 



i.e. it has a local minimum. Here m and M are two mass scales, which we will discuss in 
more detail below. The full kinetic function is sketched in figure 2. When the scalar reaches 
the minimum, the vacuum shifts, in the sense that the scalar develops a time-dependent 
vacuum expectation value that is growing linearly with time, (0) = m 2 t. The Lagrangian, 
perturbed to quadratic order around this time- dependent vev, is proportional to 



Amazingly, the perturbations in the vicinity of the minimum (P t x ~ 0) are not ghost-like, 
despite the —A term in the Lagrangian. The energy density and pressure are given by 



£smallX — P{X)smaUX ~ M 4 (— A + A 2 ), 



(41) 



C oc (2XP X x + Px)(S<p) 2 - Px(V50) 2 . 



(42) 



p = 2P X A - P, 



p = P, 



(43) 



7 In fact, it has recently been discovered in [61] that both types of theories are closely related: they are 
virtually identical when one considers time dependence alone, the main differences residing in spatial 
gradient terms. 
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P(X) 




FIG. 2: The kinetic and potential functions for the ghost condensate bounce, as used in new 
ekpyrotic models. 

so that p + p = 2PxX. Thus, at the minimum, the ghost condensate is right on the verge 
of violating the NEC, its equation of state being that of a cosmological constant p = —p. 
But if the field is pushed to even smaller X values, where P t x < 0, then it violates the 
NEC, without the appearance of ghosts, and can cause the universe to revert smoothly from 
contraction to expansion. 

For such a scenario to work, a number of conditions must be met: first, there needs to 
be a steeply rising potential to slow the scalar down after the end of the ekpyrotic phase, in 
order for the ghost condensate phase to start - see figure 2 for a sketch of the shape of the 
potential that is required. In fact, as described in [34], the potential must rise roughly as 
V oc — -^cj) for consistency. 

Secondly, even though there are no ghosts, there are still gradient instabilities when the 
NEC is violated, as is apparent from the last term in (42). However, as shown in [34, 59], 
these are harmless as long as the bounce is fast (on the order of one e-fold), as there is then 
no time for the instability to grow. But having a fast bounce requires the parameters of 
the model, in particular the mass scales m and M, to be rather finely tuned. It turns out 
that, for an ekpyrotic phase with fast-roll parameter e, lasting for iV e & e-folds and ending at 
a potential value of V m i n , the mass scales must satisfy the hierarchy 



In words, the effective field theory cut-off scale M must be at or above the scale of the 
ekpyrotic potential, with the scale m of the ghost condensate vev many orders of magnitude 
below. Such a tuning does not seem unreasonable. 



m e 




(44) 
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There is one further requirement that the ghost condensate model must satisfy, and it 
stems from a potential problem that has recently been discovered by Xue and Steinhardt, 
and which actually threatens to invalidate all non-singular models of a bounce [63]. The 
problem arises from an unlikely source, namely an apparently doubly subdominant cosmo- 
logical perturbation mode. Let us temporarily neglect the entropy perturbations that are 
thought to source the primordial density fluctuations in the present model. Then, as dis- 
cussed at the beginning of section IV, the dominant perturbations have a blue spectrum and 
ultimately provide a negligible contribution to the observed temperature fluctuations in the 
CMB. In slightly more detail, the comoving curvature perturbation ( obeys the equation of 
motion (in Fourier space and conformal time) 

a + 2-a + k\ k = o, (45) 

z 



where z = a^2e/c 2 s . c s denotes the speed of sound of the fluctuations, and can be read off 
from (42), 

c 2 = ^ (46) 

° s 2XP XX + P X l4bj 

On large scales, i.e. for small k, the equation of motion is solved by 

& = ^ + W$' (47) 

where c±, c 2 are constants of 0(1), as can be seen by matching to the Minkowski vacuum in 
the far past. The first mode is the one with the blue spectrum n s = 3, and which can be 
neglected. The second mode has an even bluer spectrum, n s = 5, and is typically not even 
mentioned. But it is instructive to re-write the integral in (47) as 

Cl nt fdr fdtc 2 f dtc 



f*L= f^= ! dtCs (48) 
~ J z 2 J 2a 3 e J 3a 3 {l + wY 1 ; 



c 2 Vk 

During the bounce, we must have w < — 1, and so it looks like the integral might blow 
up in the approach to the bounce phase, which would signal a breakdown of perturbation 
theory. However, the speed of sound also goes to zero as w reaches —1, and so the integral is 
finite. It is nevertheless dominated by the contribution close to the bounce phase, and can 
be evaluated using (43), (46) and the approximation of neglecting the friction term in the 
equation for X, X — V^y/X, 

Cg f 2XP X -P dX ^ -1 

c 2 Vk J 2a^X(2XP xx + P x )(-V (j) )^ V^y/X^e' 

Since the speed of the scalar at the bounce >/ X bounce is necessarily very small, we can see 
that in the approach to the bounce, this initially doubly subdominant term gets massively 
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amplified, and threatens to dominate over all other perturbation modes. As shown in [63], 
the ratio of this mode to the curvature perturbation ( s produced by the entropic mechanism 
is 

/-int i 

C* V ,<t> 

where denotes the number of e-folds of ekpyrosis remaining after mode k of ( mt has exited 
the horizon. For modes of observational interest ~ (20 — 30), while N e k > 60, so for many 
models of non-singular bounces this effect will be problematic. However, as we saw above, in 
the present model the bounce must be fast in order to avoid gradient instabilities, and this 
means that the potential must be steep, |V^| ~ 1/m 2 ^> 1. In fact, the consistency relation 
(44) implies that for the model discussed here, m 2 e Nek is always small, and so the doubly 
blue mode remains subdominant. Hence, the condition that the bounce be fast avoids the 
onset of all known potential instabilities. 

From a model-building point of view, it is very satisfying to see that a consistent non- 
singular bounce model can be constructed 8 . What is not so clear yet is whether such a model 
could also arise from string theory. There exist some general arguments to the contrary [65, 
66], but history has shown that no-go arguments have a tendency to be overcome. In a recent 
effort to clarify this question, the ghost condensate model has been supersymmetrized [67], 
which may help in elucidating the connection to string theory. It would also be interesting 
to see if non-singular bounce models can be constructed using the closely related galileon 
theories. These theories do not suffer from gradient instabilities, and so might lead to more 
general bounce models. The relation of these models to string theory is equally unsettled at 
present. They are thought not to arise in string theory by the arguments of [65], but on the 
other hand, they can be derived, very suggestively, from the dynamics of branes [68] . They 
have also been supersymmetrized recently [61], and so these questions may find a clarification 
soon. 



C. Nonperturbative Bounces 

The AdS/CFT correspondence is thought to provide a non-perturbative definition of 
a theory of quantum gravity. This makes one wonder whether, in this framework, one can 
finally tackle the question of cosmological singularities in a full quantum mechanical way, and 
study bounces non-perturbatively. In this vein, Craps et al. have attempted to construct 
models of a crunch of Anti-de Sitter space [69-71], which are described in detail in his 

8 See the review [64] for other attempts at constructing bounce models. 
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contribution to the present volume. Hence, we will limit ourselves to just a few remarks 
here, to point out the main differences with the colliding branes bounce discussed above: 
the AdS crunch corresponds to a high curvature crunch, and thus can be expected to lead 
to overproduction of particles at the crunch, with subsequent rapid re-collapse. Numerical 
studies in the context of a 5-dimensional model confirm this expectation [72], but there is the 
hope that a better behavior of the coupling constant in the dual conformal field theory for 
the 4-dimensional case might evade this problem. Also, in the models studied so far, there is 
no ekpyrotic matter present, and thus there is no matter component that can halt the rapid 
growth of anisotropies. Thus, the models in question may describe the quantum version of 
chaotic mixmaster collapse. This in itself will be very interesting, but at the same time it 
will also be of interest to see if the framework can be extended to low-curvature crunches 
more like that of compactified Milne space. 

D. Reheating 

In all models of a cosmic bounce, the ordinary radiation and matter content of the universe 
is thought to arise at the bounce itself. For singular bounces, one expects the brane collision 
to be slightly inelastic, so that a small part of the collision energy will be converted into 
radiation and matter [11]. Note that in that case the scalar fields do not disappear by 
decaying into radiation and matter, as is typically assumed in inflation. Rather, the scalar 
fields become cosmologically relevant again at later times, when they act as quintessence, 
and their non-decay is essential in the construction of a cyclic model of the universe, as 
detailed in the next section. For non-singular bounces on the other hand, the reheating 
mechanism most discussed so far involves the decay of the field governing the bounce, e.g. 
the ghost condensate field, into radiation and matter just after the bounce [34]. Such a 
mechanism is most suitable for single-bounce models of the universe, but appears more 
difficult to incorporate into a cyclic model. 

It seems fair to say that the topic of reheating has not been studied in sufficient detail and 
generality in the context of cosmic bounces; this is certainly an area that deserves further 
research. 

VI. COSMIC CYCLES, THE PHOENIX AND PREDICTIVITY 

By having a contracting phase before the current expanding phase, have we not simply 
pushed the question of initial conditions back to an earlier time, in much the same way as 
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FIG. 3: The effective potential in the cyclic universe, with the different phases of evolution. 

the theory of inflation does, though with the added complication of a bounce? Then all 
we have gained is a new mechanism for generating cosmological perturbations (though we 
already had one with inflation), but no answers to any other basic cosmological questions. 
However, imagine that we link up alternating contracting and expanding phases into a cyclic 
universe, could it then be that one cycle prepares the initial conditions for the next one? 
Or, assuming that there was some original big bang-like event, but which did not produce a 
universe like ours - could subsequent cycles evolve more and more into universes like ours? 
If so, such an attractor framework might be highly predictive, much more so than a theory 
which produces all types of universes, but with no clear preference for one like ours 9 . The 
cyclic model which we will discuss now is an attempt to build precisely such a model. 

The cyclic universe, first proposed by Steinhardt and Turok [11], posits that the current 
dark energy phase that we find ourselves in will be followed in the distant future by a 
contracting ekpyrotic phase, which will lead to a bounce and a new phase of radiation and 
matter domination, followed once more by dark energy, and so forth. As an effective theory, 
this can be described by a scalar field moving in the potential shown in figure 3. The figure 
shows the potential only in one direction, and one should keep in mind that other, transverse, 
directions can also be relevant: in particular, during the ekpyrotic phase, there is a transverse 
unstable potential, as described in section IV, which will play an important role here. 

The cyclic model is formulated as an effective theory, without specifying a precise mi- 
crophysical origin, although the most precise formulation to date is based on the colliding 



9 This kind of argument can be made precise, see for example [73] . 
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branes picture stemming from heterotic M-theory 10 . It is useful to briefly go over the different 
phases of the model, and to follow the evolution of the scale factor and the Hubble parameter 
[74]. Let us start with the current dark energy phase, when the scalar field takes its largest 
values. During this phase, the Hubble parameter is roughly constant and the universe grows 
by N de e-folds as the scalar field is very slowly rolling to smaller field values. Eventually 
the potential becomes negative, the field speeds up, the universe starts to contract again 
and a new ekpyrotic phase takes place. From section III we know that the scale factor now 
remains roughly constant, while the Hubble parameter grows by e Nek = V^J^/V^ 2 , where 
V denotes the current dark energy density. The ekpyrotic phase comes to an end when the 
scalar reaches the bottom of the potential at V min) and is followed by the bounce phase. If 
we assume the colliding branes model of the bounce, the dynamics will be dominated by 
the kinetic energy of the scalar, with the scalar shooting off to — oo and bouncing back after 
the brane collision. By the end of the bounce/kinetic phase, the universe has grown by the 
modest amount V^/T^ 3 , where T r denotes the temperature to which the universe reheats 
at the brane collision. Meanwhile, the Hubble parameter has shrunk by a factor T^/V^?- 
After the bounce, the scalar field flies back over the potential well, slows down due to Hub- 
ble friction and comes to a halt on the plateau at large field values. Right before the dark 
energy phase, while the energy density of the universe is still dominated by radiation and 
matter, the universe grows by a factor of about T r /T , with the Hubble parameter shrinking 
by T 2 /T r 2 , where T = V^ 4 is the current temperature of the universe. Thus, putting every- 
thing together, the Hubble parameter returns to its starting value after one cycle (as do all 
such locally measurable quantities), while the scale factor grows by the large amount of 



over the course of one cycle. This net growth is essential to the viability of the cyclic universe. 
For one, it solves a puzzle about the second law of thermodynamics, which asks how a cyclic 
universe is possible, given that the entropy must increase all the time. The resolution here is 
that, although the total entropy does increase continually, the entropy density gets diluted 
to the same small value again after each cycle. 

We should now go back to the ekpyrotic phase, and address the instability of the potential 
that is necessary for the generation of perturbations. Cyclic universe and instability are two 
concepts that do not seem to go very well together a priori However, this instability may 

10 It is at present an open question whether a consistent cyclic model of the universe can also be constructed 
using a non-singular bounce of the ghost condensate type, for example. 




(51) 
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actually lie at the heart of the high level of predictivity of the model. As described in section 
IV, during the ekpyrotic phase the background trajectory of interest is the one that rolls 
unstably along a ridge in the potential, with the potential falling off steeply both left and 
right. For slightly different initial conditions, the ekpyrotic phase will be of short duration: 
the trajectory will veer off to the side, ekpyrosis and the associated flattening of the universe 
will end prematurely, and the corresponding region of the universe will collapse in a chaotic 
mixmaster crunch. At the end of such a crunch, we will conjecture that a gravitationally 
collapsed region remains, in the form of a black hole presumably, so that this region will stop 
growing. Only a region with very special field values at the start of the ekpyrotic phase will 
make it all the way down the ridge, through the bounce and on to the next expanding phase. 
However, that region, unlike all the other collapsed regions, will grow by the enormous factor 
(51) until the start of the next ekpyrotic phase. Moreover, for N de > 60, it turns out that 
this newly expanded universe is bigger than its originator universe one cycle before [75]. 
Thus, as long as the dark energy phase lasts sufficiently long (sufficiently long meaning at 
least 600 billion years), the cyclic universe grows from cycle to cycle and selects regions of 
the universe with just the right initial conditions to produce a large, flat and habitable space 
during the next cycle. In this context, it is wrong to say that the ekpyrotic phase requires 
special initial conditions. Rather, this phoenix universe [75, 76], reminiscent in spirit of the 
phoenix universe discussed by Lemaitre [77] and Dicke, produces a vast range of "initial" 
conditions before the onset of the ekpyrotic phase 11 . Almost the entire universe will then 
collapse irredeemably, except for some small regions, which are necessarily exceptionally flat, 
homogeneous, isotropic but with small quantum generated density perturbations, growing 
into the kind of universe we see around us. And should the universe that is produced in such 
a way happen to be small (compared to the presently observed Hubble volume, say) then 
this is hardly an obstacle, as each subsequent cycle will make it bigger and bigger. Unlike 
eternal inflation, the phoenix universe renders the vast majority of space habitable. Note 
also that it is not necessary for the cycles to have occurred indefinitely into the past. A 
beginning of the universe, via a tunneling event for example, is entirely conceivable in this 
framework, as long as such an event produces one region with the right conditions for the 
ekpyrotic phase to proceed successfully. Whether the universe had a beginning or not is 
unfortunately difficult to establish experimentally in this framework. For better or worse, 
the instability implies a near total isolation from even the immediately preceding cycle, with 



Note that, given the assumed shape of the potential, this range of "initial" conditions does not however 
include conditions which would allow eternal inflation to occur. If this were the case, we would of course 
face the measurement problem once again. 
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FIG. 4: The global structure of the cyclic universe. 



only the quantum generated fluctuations surviving the bounce. Any small excess energy 
present at the onset of the last ekpyrotic phase in the region that is supposed to evolve into 
our observable universe would automatically cause the corresponding field trajectory to fall 
off the potential ridge and prevent that region from making it through the bounce. On the 
positive side, this means that there is no paradox in us not having found messages across 
the sky from far-advanced civilizations of the last cycle(s). 

As speculative as some of these discussion may sound, the cyclic universe described above 
leads to concrete pre- and postdictions, which can be confronted by experiment. The most 
important predictions, regarding non-gaussianities in the CMB and gravitational waves, were 
described in section IV. But there is also an important postdiction, concerning the amplitude 
Q of the primordial density perturbations [78] . In other models of the early universe, notably 
inflationary models, this quantity is always adjusted by hand to fit the observed value of 
Qobs ~ 1CT 5 , even though the observed value is typically not favored from a theoretical 
point of view. By contrast, the cyclic model (with a bounce modeled by colliding branes) 
dynamically selects regions of the universe with the observed value of Q\ This is entirely 
analogous in spirit to the selection of flat regions of the universe discussed in the previous 
paragraph. Since this is a significant result, it is worth deriving it in a little bit more detail. 
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The ekpyrotic phase evolves according to the scaling solution (22), or, using conformal time, 

a = {-r) l '\ ^ = V2/(V~er), V = -eU\ (52) 

with the conformal Hubble rate % = a. The subsequent kinetic phase is characterized 
by 

a={-rfl\ g = V3/(V2r). (53) 

Matching the conformal Hubble rates in the two phases at the time that the kinetic phase 
begins, implies 

-2r fc = y/e/Vmin. (54) 

Recalling also from the discussion below Eq. (40) above that the inter-brane distance is 
given by 

d = 2y \Tk\, (55) 

we get the relation V m in = ^Uo/d 2 . An important additional information is that the extension 
of the solution (40) to the ekpyrotic phase implies that the inter-brane distance is approxi- 
mately constant during that phase, with the value d given above (see [48, 78]). Hence, the 
amplitude of the curvature perturbations generated by the entropic mechanism, given 
above in (26), can be re-expressed very simply as 

Qc ~ w^d ~ 7' (56) 

where we have taken e ~ O(10 2 ) both because this is consistent with observations of the 
spectral index, and because, from a theoretical point of view, the lowest allowed values for 
e correspond to the least amount of fine-tuning in the potential. For d, we can assume 
the value d « 10 35 suggested by Hofava-Witten theory, and which implies the observed 
value of Newton's constant. The crux of the argument comes from considering y , the brane 
collision velocity. As discussed in section V, it cannot be arbitrarily high, as this would 
lead to an overproduction of matter (wrapped M2-branes) at the brane collision, with a 
subsequent rapid collapse of the corresponding region. In other words, if y is relativistic in 
some region, then that region undergoes gravitational collapse at the bounce, and becomes 
irrelevant regarding the large-scale evolution of the cyclic universe in much the same way as 
those regions with the wrong initial conditions for ekpyrosis described above. Thus, the vast 
majority of the universe must have y < 0.1, or 

Qc < ID- 4 ' 5 . (57) 

Regions in which the brane collision velocity is higher lead to a larger and thus to more 
structure, making it highly likely to find oneself in a region of the universe with an observed 
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value of Q as close as possible to the upper bound (57), see also figure 4. It is certainly 
intriguing that the observed value of Q agrees so closely with this postdicted value! 

To conclude, it seems worthwhile to dwell for a little longer on the main theme discussed 
in this last section, namely the contrast between the instability of the ekpyrotic phase and the 
net expansion over the course of each cycle. The discarding of large regions of the universe 
at the end of each cycle provides a type of dissipation, which is just what is needed to have 
attractor behavior of the model. Hence, perhaps somewhat counter-intuitively, it is precisely 
the interplay between fragility and amplification that lends the model its high predictive 
power. If it was easy to make all sorts of large, reasonably flat universes, with all sorts of 
density perturbations, all allowing for the formation of gravitationally bound structures, it 
would be very difficult to understand the broad features of our universe. We would have 
to conclude that our part of the universe is as it is because of sheer coincidence. Before 
accepting such a conclusion though, we should test those models in which the specialness of 
our universe is not due to mere historical accident, but due to mathematical and physical 
necessity. 
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